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1. Introduction 

Long ago it was conjectured by Matsui and Satz [|l|] that melting of different quarkonia states 
due to color screening can signal Quark Gluon Plasma formation in heavy ion collisions. The 
problem of quarkonia dissolution has been studied using potential models [[|^|,^], [7|]. However 
it is very unclear if such models are valid at the finite temperature [Jsj] . Therefore the problem of 
dissolution of different quarkonia states should be studied in terms of the corresponding meson 
spectral functions. For charmonium such studies appealed only recently and suggested, contrary to 



potential models, that J/y and rj c survive till temperatures as high as \ .6T C [Q, 10, 11]. It has been 



also found that % c melts at temperature of about 1.1T C [11]. It is expected that there is a hierarchy of 
dissolution temperatures for different quarkonia, states corresponding to hierarchy of their sizes, i.e. 
we expect a sequential dissolution pattern where larger states dissolve at smaller temperatures. In 
this contribution we present our preliminary study of bottomonia spectral functions and correlators 
at finite temperature. 

2. Lattice correlators and spectral functions 

In our lattice investigation we calculate correlators of point meson operators of the form ^H/, 
where T is one of the Dirac matrices which fixes the quantum number of the channel. We consider 
pseudo-scalar, vector, scalar, and axial-vector channels which correspond to rjb, T, Xbo and %b\ 
respectively. The Euclidean time correlator G{x,T) we calculate on the lattice is an analytic con- 
tinuation of the real time correlator, G(x, T) = D > (—iz, T). Due to this relation the meson spectral 
function, defined as 

D>((0,T)-D<(co,T) 
a(co,T) = , (2.1) 
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is related to the Euclidean correlator through the integral representation 

G(r,T) = J dcoo((0,T)K(T,(0,T) (2.2) 

c OS h(fl)(T-l/(2r)) 
* (T ' W ' r)= smh(fl,/(2r)) • (2 " 3) 

Thus to reconstruct the spectral function from the lattice correlator G(z, T) this integral representa- 
tion should be inverted. Since the number of data points is less than number of degrees of freedom 
(which is £?(100) for reasonable discretization of the integral ) spectral functions can be recon- 



structed only using the Maximum Entropy Method (MEM) []12p. In order to have sufficient number 
of data points either very fine isotropic lattices | fU| ] or anisotropic lattices [^, [l0|] should be used. 
Another difficulty arises due to the large quark mass, discretization errors 6{am c ^) are present 
in the heavy quark system. To remove these discretization errors we use the Fermilab approach 



We performed calculation of bottomonia correlators in quenched QCD using anisotropic lat- 
tices and Wilson gauge action. We used three value of the gauge coupling j8 = 5.9, 6.1, 6.3 and 
anisotropy t, = a s /a t = 4. When the lattice spacing is set using the Sommer scale ro = 0.5 fm the 
above gauge coupling correspond to lattice spacing = 5.88, 8.18, 10.89 GeV. Further details 



of the lattice action together with the parameters can be found in Ref. [15]. 
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Figure 1: Bottomonia spectral functions in pseudo-scalar (left) and scalar (right) channel for different lattice 
spacing. 



To be able to reconstruct the spectral functions high statistical accuracy for bottomonia corre- 
lators is required. This makes the calculations computationally intensive, so we performed them 
on a prototypes of RBRC QCDOC machine. It is a dedicated Lattice QCD machine developed 
by physicists from Columbia University, BNL, RIKEN and UKQCD. Three such machines, each 
reaching about lOTFlops peak performance, are currently installed at BNL and EPCC. Our pro- 
totypes were single-motherboard machine at about 50 GFlops peak. Such resources are still not 
adequate for the full QCD simulations, so we use quenched approximation, which is equivalent to 
neglecting quark loops. Typical statistics gathered was 500 to 1000 measurements, separated by 
400 updates. 



3. Spectral functions at zero temperature 

Bottomonia spectral functions at T = are shown in Fig. [T] for the pseudo-scalar and scalar 
channels. The first peak is independent of the lattice spacing and corresponds to T]/, and Xb respec- 
tively. The position of the higher peaks depends on the lattice spacing, so they cannot be interpreted 
as physical states. Similar results have been found for the vector and axial-vector channels. 



4. Correlators and spectral functions above deconfinement 

We would like to know what happens to different bottomonia states at temperatures above the 
deconfinement temperature T c . With increasing temperature it becomes more and more difficult to 
reconstruct the spectral functions as both the number of available data points as well as the physical 
extent of the time direction (which is l/T) decreases. Therefore it is useful to study the temper- 



ature dependence of bottomonia correlators first. From Eq. ( |2.3[ ) it is clear that the temperature 
dependence of bottomonia correlators come from two sources: the temperature dependence of the 
spectral function and temperature dependence of the integration kernel K(z,a,T). To separate 
out the trivial temperature dependence due to the integration kernel, following Ref [jllj] at each 
temperature we calculate the so-called reconstructed correlator 

G recon (z 1 T)= d(oa((O,T = 0)K(t,(o,T) (4.1) 
Jo 
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Figure 2: The ratio G(t, T)/G rec0 n('Z, T) for different temperatures at /3 = 6.3. 



Now if we assume that there is no temperature dependence in the spectral function - then the ratio of 
the original and the reconstructed correlator should be close to one, G{t, T) / G recon (z, T) ~ 1 . This 
way we can identify the cases when spectral function itself changes dramatically with temperature. 
This gives reliable information about the fate of quarkonia states above deconfmement. In Fig. || we 
show this ratio at j8 = 6.3 for several temperatures corresponding to N t = 32, 24, 20, 16 for pseudo- 
scalar and scalar channels. As one can see from the figures for pseudo-scalar channel this ratio is 
close to one and shows no temperature dependence up to temperatures T = 2.3T C . Thus we expect 
that rjb will exist in the plasma phase with almost no modification of its properties. On contrary in 
the scalar channel the correlator shows large temperature dependence and G(z,T)/G recon (z,T) is 
far from unity already at 1.15r c . Similar behavior of the correlators was found in the axial- vector 
channel. This suggest modification or eventually dissolution of %b above deconfmement. We find 
very similar temperature dependence of the correlators at larger lattice spacings corresponding to 
P = 5.9 and 6.1. This suggest that despite possible lattice artifacts in the correlators and spec- 
tral functions medium modifications of bottomonia properties at finite temperature can be reliably 
understood. 

More detailed information on different bottomonia states at finite temperature can be obtained 
by calculating spectral functions using MEM. The result of these calculation is show in Fig. ||. As 
expected the state survives in the deconfined phase with no mass shift. In the scalar channel 
the first peak becomes broader and smaller with increasing the temperatures, signaling dissolution 
of the %b state. This behavior of %b state is quite unexpected as it has similar size and binding 
energy as the r\ c and J/y states. Thus in this respect, the expected sequential dissolution pattern of 
quarkonia states does not seem to hold. 
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Figure 3: Bottomonia spectral functions in the pseudo-scalar (left) and scalar channel (right) above decon- 
finement for j3 =6.3. 
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